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Abstract 

The effect of a high-frequency microwave irradiation on the averaged sta- 
tionary charge in a single-electron tunneling box (SET-box) with super- 
conducting island is investigated. The resulting {Q)~Qg characteristics are 
modified compared with those from the well-known constant voltage-driven 
box. It can be shown that there is a cross-over from the 2e-periodic regime 
in Qg corresponding to the parity effect to le-periodic behavior not only 
with increasing temperature but also with increasing irradiation amplitude. 
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I. INTRODUCTION 



The tunneling of single electrons (SET) through ultrasmall junctions has become of 
much experimental and theoretical interest (see e.g. ^)). The main conditions to have 
charging effects (Coulomb blockade) are low temperatures (fc^T <^ Ec) and large tunnel 
resistances {Rt ^ Rq = h/e^ = 25.8 kVL). This guarantees that single-electron effects 
are lifted out of thermal and quantum fluctuations, e and h are the elementary charge 
and Planck's constant, respectively, /c^ is Boltzmann's constant. Ec denotes the relevant 
Coulomb energy under consideration. 

We consider the so called SET-box formed by a SET junction and a capacity Cs in 
series as described in Fig. 1. It is known that the mean electric charge on the island of 
this SET-box can be increased in integer steps by a continuous increase of the external 
voltage V = Vo (see e.g. 0) corresponding to an increase of the charge Qg at the capacity 
Cs {Qg = CsVo)- The staircase function {Q) vs. Qg is the sharper the temperature 
is lower (see Fig. 2). The tunnel junction is characterized by the capacity C and the 
imaginary part of the quasi-particle current amplitude Imlg which is related with the 
phenomenological tunnel resistance Rt- 

A new aspect which has thoroughly been investigated as well theoretically and exper- 
imentally are parity effects on the superconducting island [3-11]. It has been observed 
that the parity based 2e-periodicity of the charge characteristics disappears in the vicin- 
ity of a cross-over temperature which is much lower than the critical temperature of the 
superconductor. 

There are some articles dealing with time-dependent phenomena in mesoscopic tun- 



nel junctions [12-16]. We have investigated in an earlier paper [|17| the influence of a 
microwave irradiation on the stationary behavior of a normalconducting double-junction. 

The aim of this paper is to investigate the influence of a microwave irradiation on 
SET-box parity effects in stationary regime. Because of the fact that in some sense the 



microwave irradiation acts similarly like temperature increase the cross-over mentioned 
above can also be achieved by a sufficient strong irradiation. 



II. NORMALCONDUCTING BOX 

For a normalconducting island the relevant energy function which enters the transition 
rates reads as 

n labels the number of electrons on the island. The index "ch" means that this energy 
is the charging energy. The charge Qg at the capacity Cs is connected with the external 
voltage Vo via the relation 

The master equation describing the dynamics of the occupation probabilities of the 



SET-box reads as 0.19 



^(^nl^) — ~(^n-l,n + ^ n+l,n)o'ni't) + r„_„+iCr„_|_i (t) + T n,n-l'^n-lit) ■ (2) 

The rates Tn±i,n which depend on the energy differences (cf. Eq. (|I|), Ec = e^/ {2{C+Cs))) 

^KliAQs) = KiiiQ,) - K\Qa) = Ec{l±2nT 2QJe) , (3) 
are given by 

TniiAQa) = llmI,[-AE^,l,^^{Qg)/h] {l - eM^Kii,niQ9)/kBT]y' . (4) 

Iq are the quasi particle current amplitudes of the NIN junction. The rates r„j-i „ cor- 
respond to the tunneling processes connected with the charge transitions n nil, 
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respectively. The succession of indices of the rates is the same as for the energy differ- 
ences, the first index refers to the final state and the second to the initial state. The mean 
stationary charge in the box is given by 

(Q) = e^w„(Qg) , (5) 



where cx^ are the stationary solution of Eq. (^. They can be written as pO| , p! 



Jt — 1 oo 

0"n(<5g) = T7 II ^m+l,m{Qg) H ^m'-l,m'{Qg) ■ (6) 
m=-ao m'=n+l 

Af guarantees the correct normalization. It can be shown that the mean charge (|^) is 
given by a Boltzmann distribution B] 



-E^'^{Qg)/kBT 



n 

This Boltzmann distribution shows that the charge in the box corresponds to a stationary 
state and dynamical properties expressed by the current amplitudes do not play any role. 



III. NORMALCONDUCTION BOX WITH IRRADIATION 

The influence of a microwave irradiation on single-electron tunneling through single- 
junctions has been studied in Refs. [13-16]. The microwave irradiation is described by an 
additional oscillating voltage part 

V{t) = Vo + Vi cos ujt . (8) 

In this way the microwave field is treated classically neglecting possible back reactions 
due to tunneling. This is reasonable because all calculations are only done in first order 
perturbation theory. The effect of a microwave irradiation in case of the box is very 
similar to that of the double-junction because there is the island which can only be 
charged in integer quantum units. 
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But now the time averaged charge does not more satisfy a Boltzmann distribution. It 
is rather given by formula (|^) where the transition rates have been modified (cf. Eq. (^)) 
r — > r. The new rates read as 

r„±i,„(Q,) = ; E mimI,[-AE:i,JQ,)/n - ru] {l - e(^^S'i.^(««)+^'M/%T}-^ . (9) 
e J, 

The parameter a is given by 

huj 

and Ji are the Bessel functions of the first kind. The description via a Boltzmann distri- 
bution fails because the irradiation does not allow a pure stationary description. Also the 
stationary parts of the occupation probabilities carry a relic of the dynamics. From the 
mathematical point of view the reason is the sum structure of the rates which does not 
allow a factorization. This approach is based on the assumption that possible measure- 
ments of the island charge (e.g. by a on-chip SET-electrometer) would only detect the 
time-averaged charge. Therefore, these rates correspond already to an averaging with re- 
spect to the oscillations. But nevertheless they contain information from the irradiation. 
For vanishing AC amplitudes these expressions lead to the ordinary box transition rates 
without irradiation (^. 

The result is (see Fig. 3) that the increase of the AC voltage amplitude gives rise 
to a smoothing of the mean charge function compared with that for Vi = 0. In this 
way the oscillating driving voltage smears out the hard staircase and acts similarly like a 
temperature increase. This can be demonstrated in the following manner. For T = the 
rates without irradiation are nonvanishing only for AE^'^ < 0. But T > there are also 
transitions if AE'^^ > 0. In case of irradiation and T = the condition for nonvanishing 
rates is AE^^ + rhu < which means that there are also transitions for AE'^^ > 0. The 
comparison of Figs. 2 and 3 shows that there is nearly no difference for a specific choice of 
parameters. The mechanism can be understood as photon assisted tunneling via excited 
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states. The main effect is due to 1-photon processes (r = ±1) because multiple photon 
processes are strongly damped. But the specific shape of the charge function depends on 
a complicated interplay of the parameters a, uj and T. 



IV. SUPERCONDUCTING BOX 

The phenomenological treatment of the superconducting case starts with the observa- 
tion that the energy content of the island depends on the parity P„ of the electron number 
n on it. This means that for odd n there is the additional energy gap A in the energy 
function otherwise there is no gap for even n , P„ = n mod 2 . This corresponds to the 
substitution 

(g,) ^ p„(g,) = E^:^\Q,) + A ■ p„ . (10) 

The same approach as in Sec. can be done using the energy function (|TUI). This 
consideration is correct for zero temperature. Then the charge in the box jumps to the 
next integer value at the crossing point of two neighbouring En{Qg) parabolas (see Fig. 4). 
The shape of the resulting charge staircase function depends on the relation A between 
the energy gap A and the Coulomb energy Ec (A = A/Ec) (see Fig. 5). The {Q)-Qg 
characteristics shows for A > a 2e-periodic behavior in Qg. But this approach cannot be 
correct for T > because it does not explain the cross-over from the 2e-periodicity to le- 
periodicity at a cross-over temperature T* known from experiments [Q. This cross-over 
arises due to a temperature dependent difference of the free energies 6F{T) of the island 
in the even and odd state SF is reduced from A at T = by entropy contributions. 



It is not an effect of the temperature dependence of A itself because this temperature is 
still much lower than the critical temperature Tc- Therefore A = A(0) throughout this 
paper. A first order approximation for 6F{T) reads as 

5F{T) ^ A - knT In N^jf . (11) 



Neff is the effective number of states avaiable for excitation. It depends only weakly on 
temperature 

iVe// = 2iV(„)(0)VAe^/'=^^iri(A/A:BT) . 

^(n)(0) is the normal density of states at the Fermi level and V corresponds to the island 
volume. Ki labels the modified Bessel function. The cross-over temperature is defined 
by 

5F{T*) = . 

For their specific sample Tuominen et al. found 

SF{T) ^ A(l - ^) 

with T* ^ 300mi^ much less than Tc- Of course this is reasonable only for < T < T* , 
the function 6F{T) is positive definite. For < Qg/e < (1 + A)/2 and ksT < Eq one 
can use the approximation that the system is mainly governed by the two states n = 
and n = 1 (cf. Fig. 4). Generalizing the approach of Sec. (|Tp the stationary occupation 
probabilities obey the Boltzmann distribution 



a(0) 



^_(^Ef+5F{T))/kBT 



^^^^ e-Et/ksT _^ ^-{Ef+5F{T))/kBT ' ^^'^^ 



Then (Q) is given by 



(g) =e^_^,,^^^^^^^y^^_^^,^^^^^^^^^^^^^^^^ . (13) 

In Fig. 6 the box charge has been calculated over the main Q^-interval for a sample 
of temperatures using Eq. ([T3|) . There are two mechanisms how temperature changes 



the shape of the charge characteristics. The kBT-ieims, in Eq. (p!3D make the curve 



smoother but they do not shift the critical values of Qg where the mean charge reach half 
integer values. This is just reached by the 5F(T)-term. At A/ksT = 8 corresponding 
approximately to the cross-over temperature the 2e-periodicity has already disappeared 
leaving only the le-periodicity. The property {Q){Qg = 1/2) = 1/2 means that there are 
no more even-odd-effects. 

But there is another approach to that problem discussing in more detail the transition 
rates |]10|- The transition n = to n = 1 is governed by the rate 



ri,o(Q,) = ^Im/,J-Ai?fo(Q.)/^] {l - e^^?o(«.)A.T|-i ^^^^ 

This corresponds to Fi o in Eq. for a NTS junction. The current amplitude Iq^ corre- 
sponds to the NIS junction and brings the gap energy A up. But the back rate describing 
the transition n = 1 to ri = consists of two terms 

ro,i = 7o,i + rS, (15) 

where the second rate on the right corresponds again to Fq,! in Eq. (^) and reads as 

riim = llmIqA-AE^MQ,)/h] {l - e^^^oMQ.y'^BT^ . (16) 

The rate 70,1 describes back tunneling of just the one (odd) electron which does not found 
any partner for pairing from the island and is proportional to N~j^jr [jlO|. This is the 
difference to the approach of Sec. 2. Of course 70,1 depends on the energy differences 

AEQ i{Qg) too. The rates are given by "golden rule " integrals |T0|1 

00 

ri,o = y / rf6iV(„)iV(,)(e - A£;^J)|T|V(e)(l - /(e - AE'^'},)), (17) 

00 

70,1 = fj rfeiVHiV(,)(6 + AEo1)|Tp(l - f{e))f{e + AE^,^)). (18) 



A''(„) and N(^s) are the densities of states of normalconductor and superconductor, respec- 
tively. T is the tunneling amplitude. The distribution / is defined by 
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deiV(,)(e)/(e) = 1 



corresponding to the assumption that there is just one unpaired electron on the island. 
Now the mean stationary charge reads as (0 < Qg/e < 1) 



■7o,i + rg(Q,) + ri,o(Q, 



(Q) = e — • (19) 



What is the connection between Eq. (|T3D and Eq. (|19])? Using Eq. (0), Eq. (|T3D can 
be transformed in the following form 

+ e 

This expression is a rather good approximation of Eq. (|19D because for < Qg/e < 
(1 + A)/2 there is 

ri,o(Qg) oc exp(-(AEfo(Q,) + A)/A;bT)) 

and Fq"] <^ 7o,i. Also for (1 + A)/2 < Qg/e < 1 this formula yields the correct charge 
value {{Q) /e ^ 1). The explicit consideration of the possibility that just the one unpaired 
electron can tunnel back corresponds to the free energy argument in the other approach. 

V. SUPERCONDUCTING BOX WITH IRRADIATION 

The discussion of the transition rates is necessary because the treatment of irradiation 
is not possible in the Boltzmann distribution approach. But the rates ([1^, |I5|) have to be 
generalized with respect to arbitrary n and Qg. 

Tnil.n = 7n±l,n ' Pn + r^2l,n ' (21) 

Ti%AQ,) = ;Im/,[-AEf^,jg,)/^] {l - e'^^'rH.JQ.y'^sTy^ . (22) 



The rate 7 in Eq. ( pT]) describes either back tunneling of just the one (odd) electron (see 
Sec. 4) or the direct transition of just the pairing partner of the odd electron on the island. 



The other rates T describe transitions connected with the production of excited states on 
the island. 

The treatment of the case with microwave irradiation is analogous to that of Sec. 3. 
The charge is given by formula where the occupation probabilities a are determined by 
modified transition rates T. These modified rates have been phenomenological constructed 
analogous to Eq. (||)) and read as 

f ...(Qs) = E JHa)T„XAE':'{Q,) + rhoo) , (23) 



where the symbol F.,, refers to the rates (|2lD . Now the occupation probabilities cr(n) have 



to be calculated. This can be done using standard techniques (see e.g. Ref. Note 



that J2n'^{^) = 1- A sample of charge functions with different AC amplitudes has been 
plotted in Fig. 7 for fixed temperature. 



VI. DISCUSSION 

The investigation shows that the photon assisted tunneling via excited states leads 
to an AC amplitude dependent cross-over from the 2e-periodic regime to le-periodicity. 
The main effect is again due to 1-photon processes (r = ±1). The critical value of Qg 
where the mean charge reachs the value 1/2, shows that for fixed temperature T < T* (T* 
is the cross-over temperature in case on no irradiation) this value goes with increasing 
microwave amplitude to 1/2. This means that the cross-over from 2e-periodic behavior to 
le-periodicity can be reached not only by increasing temperature but also by increasing 
irradiation amplitude. An analytical formula for the cross-over amplitude cannot be 
given. The numerical simulation shows that for < j"* the le-periodic behavior is 
already established approximately for a ~ 0.1. By increasing the AC amplitude it can 
also be recognized (cf. Fig. 7) that the mechanism of "smoothing" the staircase function 
is much more complicated than that of finite temperature. The graphs for medium AC 
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amplitudes show kinks at several values Qg which are relics of the kink of the NTS current 
amplitude at the gap energy. 

The master equation approach demands that the rates should be much less than 
the other relevant frequencies. With respect to the Coulomb frequency Ec/h this is in 
principle guaranteed by Rq ^ Rt- But this should also be ensured with respect to u. 
Therefore, we have assumed that u is approximately of the same order of magnitude as 
the Coulomb frequency. If the Coulomb energy corresponds approximately to 1 K the 
radiation frequency should be in the range of 20 GHz. It is reasonable to choose u in such 
a way because this yields the natural physical frequency scale and guarantees more clearly 
effects. Up to now there are no experimental measurements known. On the contrary one 



tries to suppress photon induced processes by filtering p3|| . 
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FIGURES 



FIG. 1. Scheme of the considered SET-box 

FIG. 2. Charge in the NC box for different temperatures Ec/ksT = 100, 10, 5, 1. The 
nearly straight hne belongs to Ec/ksT = 1. 

FIG. 3. Charge in the NC box with irradiation at fixed temperature Ec/kBT = 100 for 
different AC amplitudes a = eVi/htJ = 0, 0.1, 0.5, 1, (hio/Ec = 1). 

FIG. 4. The En{Vo) parabolas for A/Eq = 0.8. The numbers n of the respective parabolas 
are given by the x-coordinate of their apices. The lowest crossing points determine the values 
where charge jumping occurs for T = (cf. Fig. 5). Note, that for A/Ec > 1 the states with 
odd number n do not participate; only two-electron tunneling processes (Andreev reflection) 
take place. 

FIG. 5. Charge in the SC box corresponding to the case of Fig. 4 for zero temperature. 

FIG. 6. Charge in the SC box in case of A/Ec = 0.8 at different temperatures 
A/ksT = 8, 10, 20, 50, 100, oo, (A^(O)VA « 10"^ corresponding to the sample of Ref. [4]). 

FIG. 7. Charge in the SC box with irradiation in case of A/Eq = 0.8 at fixed temperature 
A/ksT = 100 for different AC amplitudes a = eVi/hO = 0, 0.01, 0.02, 0.05, 0.1, {hu/Ec = 1 
, A^(0)VA ^ 10~^ corresponding to the sample of Ref. [4]). 
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